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A black hole accretion may have both the Keplerian and the sub-Keplerian com- 
ponents. The Keplerian component supplies low-energy (soft) photons while the sub- 
Keplerian component supplies hot electrons which exchange their energy with the soft 
photons through Comptonization or inverse Comptonization processes. In the sub- 
Keplerian flow, a shock is generally produced due to the centrifugal force. The post-shock 
region is known as the CENtrifugal pressure supported BOundary Layer or CENBOL. 
We compute the effects of the thermal and the bulk motion Comptonization on the soft 
photons emitted from a Keplerian disk by the CENBOL, the pre-shock sub-Keplerian 
disk and the outflowing jet. We study the emerging spectrum when both the converg- 
ing inflow and the diverging outflow (generated from the CENBOL) are simultaneously 
present. From the strength of the shock, we calculate the percentage of matter being 
carried away by the outflow and determined how the emerging spectrum depends on the 
the outflow rate. The pre-shock sub-Keplerian flow was also found to Comptonize the 
soft photons significantly. The interplay among the up-scattering and down-scattering 
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effects determines the effective shape of the emerging spectrum. By simulating several 
cases with various inflow parameters, we conclude that whether the pre-shock flow, or 
the post-shock CENBOL or the emerging jet is dominant in shaping the emerging spec- 
trum, strongly depends on the geometry of the flow and the strength of the shock in the 
sub-Kcplcrian flow. 

Keywords: accretion disk, black hole physics, shock waves, radiative processes, Monte- 
Carlo simulations 

1. Introduction 

It is well known (Chakrabarti 1990, hereafter C90) that the flow velocity is the 
same as the velocity of light c as the matter enters through the event horizon. 
However, the sound speed is never so high. Thus the incoming flow on a black 
hole is always supersonic and thus these solutions are likely to be most relevant 
in the study of the physical processes around black holes. As the flow begins its 
journey sub-sonically very far away, and becomes supersonic on the horizon, the 
flow is also known as a transonic flow. In the context of the spherical flows, Bondi 
(1952) solution of accretion and Parker (1959) solution of winds are clear examples 
of transonic flows. But they have only one sonic points. In presence of angular 
momenta, the flow may have two saddle type sonic points with a shock in between 
(C90, Chakrabarti, 1996). The solutions with shocks have been extensively studied 
in both the accretion and the winds even when rotation, heating, cooling etc. are 
included (Chakrabarti, 1990, 1996). The study demonstrates that the accretion and 
the winds are inter-related - the outflows are generated from the post-shock region. 
Subsequently, in Chakrabarti (1999, hereafter C99), Das & Chakrabarti (1999) and 
Das et al. (2001), the mass outflow rate was computed as a function of the shock 
strength and other flow parameters. Meanwhile, in the so-called two component 
advective flow (TCAF) model of Chakrabarti & Titarchuk (1995) and Chakrabarti 
(1997), the spectral states were shown to depend on the location and strength of 
the shock. Thus, C99 for the first time, brought out the relationship between the 
jets and outflows with the presence or absence of shocks, and therefore with the 
spectral states of a black hole candidates. This paves the way to study the relative 
importance between the Compton cloud and the outflow as far as emerging spectrum 
is concerned. 

Computation of the spectral characteristics have so far concentrated only on the 
advective accretion flows (Chakrabarti & Titarchuk, 1995; Chakrabarti & Mandal, 
2006) and the outflow or the base of the jet was not included. In the Monte-Carlo 
simulations of Laurent & Titarchuk (2007) outflows in isolation were used, but not in 
conjunction with inflows. In Ghosh, Chakrabarti & Laurent (2009, hereafter Paper 
I), the results of Monte-Carlo simulations in a setup similar to that of Chakrabarti 
& Titarchuk (1995) was presented. In the present paper, we improve this and obtain 
the outgoing spectrum in presence of both inflows and outflows. We also include 
a Keplerian disk inside an advective flow which is the source of soft photons. We 
show how the spectrum depends on the flow parameters of the inflow, such as 
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the accretion rates of the two components and the shock strength. The post-shock 
region being denser and hotter, it behaves like the so-called 'Compton cloud' in 
the classical model of Sunyaev and Titarchuk (1980). This region is known as the 
CENtrifugal pressure supported BOundary Layer or CENBOL. Since the shock 
location and its strength depends on the inflow parameters, the variation of the size 
of the Compton cloud, and then the basic Comptonized component of the spectrum 
is thus a function of the basic parameters of the flow, such as the specific energy, 
the accretion rate and the specific angular momentum. Since the intensity of soft 
photons determines the Compton cloud temperature, the result depends on the 
accretion rate of the Keplerian component also. In our result, we see the effects of 
the bulk motion Comptonization (Chakrabarti & Titarchuk, 1995) because of which 
even a cooler CENBOL produces a harder spectrum. At the same time, the effect of 
down-scattering due the outflowing electrons is also seen, because of which even a 
hotter CENBOL causes the disk-jet system to emit lesser energetic photons. Thus, 
the net spectrum is a combination of all these effects. 

In the next section, we discuss the geometry of the soft photon source and the 
Compton cloud in our Monte-Carlo simulations. In §3, we present the variation 
of the thermodynamic quantities and other vital parameters inside the Keplerian 
disk and the Compton cloud which are required for the Monte-Carlo simulations. 
In §4, we describe the simulation procedure and in §5, we present the results of our 
simulations. Finally in §6, we make concluding remarks. 



2. Geometry of the electron cloud and the soft photon source 

The problem at hand is very complex and thus we need to simplify the geometry 
of the inflow-outflow configuration without sacrificing the salient features. In Fig. 
1, we present a cartoon diagram of our simulation set up. The components of the 
hot electron clouds, namely, the CENBOL, the outflow and the sub-Keplerian flow, 
intercept the soft photons emerging out of the Keplerian disk and reprocess them 
via inverse Compton scattering. An injected photon may undergo a single, multiple 
or no scattering at all with the hot electrons in between its emergence from the 
Keplerian disk and its detection by the telescope at a large distance. The photons 
which enter the black holes are absorbed. The CENBOL, though toroidal in nature, 
is chosen to be of spherical shape for simplicity. The sub-Keplerian inflow in the pre- 
shock region is assumed to be of wedge shape of a constant angle <£. The outflow, 
which emerges from the CENBOL in this picture is also assumed to be of constant 
conical angle $. In reality, inflow and outflow both could have somewhat different 
shapes, depending on the balance of the force components. However, the final result 
is not expected to be sensitive to such assumptions. 
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Fig. 1. A cartoon diagram of the geometry of our Monte-Carlo simulations presented in this 
paper. The spherical inflowing post-shock region (CENBOL) surrounds the black hole and it is 
surrounded by the Keplerian disk on the equatorial plane and a sub-Kcplcrian halo above and 
below. A diverging conical outflow is formed from the CENBOL. Typical path of a photon is 
shown by zig-zag paths. 



2.1. Distribution of temperature and density inside the Compton 
cloud 

We assume the black hole to be non-rotating and we use the pseudo-Newtonian po- 
tential (Paczyhski & Wiita, 1980) to describe the geometry around a black hole. This 
potential is — 2 ( r -i) (Here, r is in the unit of Schwarzschild radius r g — 2GM/c 2 ). 
Velocities and angular momenta are measured are in units of c, the velocity of 
light and r g c respectively. For simplicity, we chose the Bondi accretion solution in 
pseudo-Newtonian geometry to describe both the accretion and winds. The equa- 
tion of motion of the sub-Keplerian matter around the black hole in the steady state 
is assumed to be given by, 

du IdP 1 
dr p dr 2(r — l) z 

Integrating this equation, we get the expression of the conserved specific energy as, 



u2 2 1 

Here P is the thermal pressure and a is the adiabatic sound speed, given by a = 
\Z"fP/p, 7 being the adiabatic index and is equal to | in our case. The conserved 
mass flux equation, as obtained from the continuity equation, is given by 

M = flpur 2 , (2) 

where, p is the density of the matter and f2 is the solid angle subtended by the flow. 
For an inflowing matter, f2 is given by, 
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where, <I> is the half-angle of the conical inflow. For the outgoing matter, the solid 
angle is given by, 

ttout = 4tt(1 - cos$), 

where $ is the half-angle of the conical outflow. From Eqn. 2, we get 

fi = a 2n ur 2 . (3) 

The quantity fi = " 7 ^" is the Chakrabarti rate (Chakrabarti, 1989, C90, 1996) 
which includes the entropy, K being the constant measuring the entropy of the 
flow, and n = is called the polytropic index. We take derivative of equations 

da 
dr 



(1) and (3) with respect to r and eliminating 4^ from both the equations, we get 



the gradient of the velocity as, 

7 1 _ 2af. 

du _ 2 (r-l)" 

dr o^-u ■ (> 

u 

From this, we obtain the Bondi accretion and wind solutions in the usual manner 
(C90). Solving these equations we obtain, u, a and finally the temperature profile 

2 

of the electron cloud (T e ) using T e = A '° fc " p , where /i = 0.5 is the mean molecular 
weight, m p is the proton mass and ks is the Boltzmann constant. Using Eq. (2), we 
calculate the mass density p, and hence, the number density variation of electrons 
inside the Compton cloud. We ignore the electron-positron pair formation inside 
the cloud. 

The flow is supersonic in the pre-shock region and sub-sonic in the post-shock 
(CENBOL) region. We chose this surface at a location where the pre-shock Mach 
number M = 2. This location depends on the specific energy e (C90). In our sim- 
ulation, we have chosen e = 0.015 so that we get R s = 10. We simulated a total 
of six cases. For Cases l(a-c), we chose rfih = 1, nid = 0.01 and for Cases 2(a-c), 
the values are listed in Table 2. The velocity variation of the sub-Kcplerian flow is 
the inflowing Bondi solution (pre-sonic point). The density and the temperature of 
this flow have been calculated according to the above mentioned formulas. Inside 
the CENBOL, both the Keplerian and the sub-Keplerian components arc mixed 
together. The velocity variation of the matter inside the CENBOL is assumed to 
be the same as the Bondi accretion flow solution reduced by the compression ratio 
due to the shock. The compression ratio (i.e., the ratio between the post-shock and 
pre-shock densities) R is also used to compute the density and the temperature 
profile has been calculated accordingly. When the outflow is adiabatic, the ratio of 
the outflow to the inflow rate is (Das et al. 2001) given by, 



D Qout ( fo \ R 



4{8(R-1) x 



3 V R 2 



3/2 



(5) 



here, we have used n = 3 for a relativistic flow. From this, and the velocity varia- 
tion obtained from the outflow branch of Bondi solution, we compute the density 
variation inside the jet. In our simulation, we have used $ = 58° and * = 32°. Fig. 



June 22, 2010 0:9 WSPC/INSTRUCTION FILE hgsgskc P rV3 



6 Ghosh, Garain, Chakrabarti and Laurent 

2 shows the variation of the percentage of matter in the outflow for these particular 
parameters. 




2 4 6 

Compression Ratio (R) 

Fig. 2. Ratio of the outflow and the inflow rates as a function of the compression ratio R of the 
inflow when the outflow is adiabatic. In our simulations, we have used the jet angle to be 58°. 



2.2. Keplerian disk 

The soft photons are produced from a Keplerian disk whose inner edge coincides 
with CENBOL surface, while the outer edge is located at 500r g . The source of soft 
photons have a multi-color blackbody spectrum coming from a standard (Shakura 
& Sunyaev, 1973, hereafter SS73) disk. We assume the disk to be optically thick 
and the opacity due to free-free absorption is more important than the opacity due 
to scattering. The emission is black body type with the local surface temperature 
(SS73): 

1/4 



T(r) » 5 x 10 7 (M^r 1/2 (M dl7 ) 1/4 (2r)- 3 / 4 



1 



K, 



(6) 



The total number of photons emitted from the disk surface is obtained by integrating 
over all frequencies (v) and is given by, 



n 7 (r) 



1^ 

he 



x 1.202057 



(T(r)Y 



(7) 



The disk between radius r to r + Sr injects dN(r) number of soft photons. 

dN(r) = 2irr6rH(r)n 7 (r), (8) 
where, H (r) is the half height of the disk given by: 



H(r) - 10 5 M dl7 



cm. 



(9) 
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The soft photons are generated isotropically between the inner and outer edge of 
the Keplerian disk but their positions are randomized using the above distribution 
function (Eq. 8) of black body temperature T(r). All the results of the simulations 
presented here have used the number of injected photons to be 6.4 x 10 s . In the 
above equations, the mass of the black hole M^h is measured in units of the mass 
of the Sun (M ), the disk accretion rate Mdn is m units of 10 17 gm/s. We chose 
M bh = 10 and Sr = 0.5r g . 

2.3. Simulation Procedure 

In a given simulation, we assume a given Keplerian rate and a given sub-Keplerian 
halo rate. The specific energy of the halo provides hydrodynamic properties (such 
as number density of the electrons and the velocity variation) and the thermal 
properties of matter. Since we chose the Paczynski-Wiita (1980) potential, the radial 
velocity is not exactly unity at r = 1, the horizon, but it becomes unity just outside. 
In order not to over estimate the effects of bulk motion Comptonization which is 
due to the momentum transfer of the moving electrons to the horizon, we shift 
the horizon just outsize r = 1 where the velocity is unity. The shock location of 
the CENBOL is chosen where the Mach number M = 2 for simplicity and the 
compression ratio at the shock is assumed to be a free parameter. These simplifying 
assumptions are not expected to affect our conclusions. Photons are generated from 
the Keplerian disk according to the prescription in SS73 as mentioned before and 
are injected into the sub-Keplerian halo, the CENBOL and the outflowing jet. 

In a simulation, we randomly generated a soft photon out of the Keplerian disk. 
The energy of the soft photon at radiation temperature T(r) are calculated using 
the Planck's distribution formula, where the number density of the photons (n 7 (S)) 
having an energy E is expressed by 



where b — l/kT(r); ((3) = ^ 3 = 1-202 is the Riemann zeta function. 

Using another set of random numbers we obtained the direction of the injected 
photons and with yet another random number we obtained a target optical depth t c 
at which the scattering takes place. The photon was followed within the CENBOL 
till the optical depth (r) reached r c . The increase in optical depth (dr) during its 
traveling of a path of length dl inside the electron cloud is given by: dr — p n adl, 
where p n is the electron number density. 

The total scattering cross section a is given by Klein-Nishina formula: 



n 7 (£) = 



1 



b 3 E 2 (e bE - 1) 



-l 



(10) 



2C(3) 




(11) 



where, x is given by, 




(12) 
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r e = e 2 /mc 2 is the classical electron radius and m is the mass of the electron. 

We have assumed here that a photon of energy E and momentum — O is 
scattered by an electron of energy 7771c 2 and momentum jJ* = 7771"^, with 

7 = [1-^1 and /_t = fJ.v. At this point a scattering is allowed to take 
place. The photon selects an electron and the energy exchange is computed through 
Compton or inverse Compton scattering formula. The electrons are assumed to 
obey relativistic Maxwell distribution inside the CENBOL. The number dN(p) of 
Maxwellian electrons having momentum between p to p + dp is expressed by, 

(13) 



dN(p) = exp{-(p 2 c 2 + n^^f^/kT^dp. 



Generally, the same procedure as in Paper I was used, except that we are now 
focusing on those photons also photons which were scattered at least once by the 
outflow. We are especially choosing the cases when the jet could play a major role 
in shaping the spectrum. 



3. Results and Discussions 

In Fig. 3(a-c) we present the velocity, electron number density and temperature 
variations as a function of the radial distance from the black hole for specific energy 
e = 0.015. rrid — 0.01 and nih = 1 were chosen. Three cases were run by varying 
the compression ratio R. These are given in Col. 2 of Table 1. The corresponding 
percentage of matter going in the outflow is also given in Col. 2. In the left panel, 
the bulk velocity variation is shown. The solid, dotted and dashed curves are the 
same for R = 2 (Case la), 4 (Case lb) and 6 (Case lc) respectively. The same line 
style is used in other panels. The velocity variation within the jet does not change 
with R, but the density (in the unit of cm~ 3 ) does (middle panel). The doubledot- 
dashed line gives the velocity variation of the matter within the jet for all the 
above cases. The arrows show the direction of the bulk velocity (radial direction in 
accretion, vertical direction in jets). The last panel gives the temperature (in keV) 
of the electron cloud in the CENBOL, jet, sub-Keplerian and Keplerian disk. Big 
dash-dotted line gives the temperature profile inside the Keplerian disk. 

In Figs. 4(a-c), we show the velocity (left), number density of electrons (middle) 
and temperature (right) profiles of Cases 2(a-c) as described in Table 2. Here we 
have fixed nid = 1-5 and nih is varied: nih = 0.5 (solid), 1 (dotted) and 1.5 
(dashed). No jet is present in this case (R = 1). To study the effects of bulk motion 
Comptonization, the temperature of the electron cloud has been kept low for these 
cases. The temperature profile in the different cases has been chosen according to 
the Fig. 3b of CT95. The temperature profile of the Keplerian disk for the above 
cases has been marked as 'Disk' 



Tabic 1 




R, P m 


"int 


JVcs 




Njct 




Neap 


P 




la 


2. 58 


2.7E+08 


4.03E+08 


1.35E+07 


7. 18E | 1)7 


8.39E | US 


3.35E + 05 


63 


0.43 


lb 


4, 97 


2.7E+08 


4.14E+08 


2.39E-- 06 


1.28E + 08 


8.58E | 08 


3.27E + 05 


65 


1.05 


lc 


6. 37 


2.7E+08 


3.98E + 08 


5.35E + 07 


4.75E + 07 


8.26E+08 


3.07E+05 


62 


-0.4 



In Table 1, we summarize the details of all the Cases results of which were de- 
picted in Fig. 3(a-c). In Col. 1, various Cases are marked. In Col. 2, the compression 
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Fig. 3. (a-c): Velocity (left), density (middle) and the temperature (right) profiles of Cases l(a-c) 
as described in Table 1 are shown with solid (R = 2), dotted (R = 4) and dashed (R = 6) curves, 
nid = 0.01 and rri/, = 1 were used. 




Fig. 4. (a-c): Velocity (left), density (middle) and the temperature (right) profiles of Cases 2(a- 
c) as described in Table 2 are shown with solid (rrih = 0.5), dotted (1) and dashed (1.5) curves, 
nid = 1.5 was used throughout. Velocities are the same for all the disk accretion rates. 



ratio (R) and percentage P m of the total matter that is going out as outflow (see, 
Fig. 2) are listed. In Col. 3, we show the total number of photons (out of the total 
injection of 6.4 x 10 8 ) intercepted by the CENBOL and jet (N int ) combined. Column 
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4 gives the number of photons (N cs ) that have suffered Compton scattering inside 
the flow. Columns 5, 6 and 7 show the number of scatterings which took place in 
the CENBOL (N cen b i), in the jet (Nj et ) and in the pre-shock sub-Keplerian halo 
(Nsubkep) respectively. A comparison of them will give the relative importance of 
these three sub-components of the sub-Keplerian disk. The number of photons cap- 
tured (N cap ) by the black hole is given in Col. 8. In Col. 9, we give the percentage 
p of the total injected photons that have suffered scattering through CENBOL and 
the jet. In Col. 10, we present the energy spectral index a (1(E) ~ E~ a ) obtained 
from our simulations. 



Tabic 2 






™d 


«int 


N cs 


jV ms 


"subfcep 


Neap 


V 


«li a 2 


2a 


0.5 


1.5 


1.08E - 06 


2.13E+08 


7.41E + 05 


3.13E | 1)8 


1.66E + 05 


33.34 


-0.09, 0.4 


2b 


1.0 


1.5 


1.22E+06 


3.37E- 118 


1.01E + 06 


6.82E | 1)8 


2.03E + 05 


52.72 


-0.13, 0.75 


2c 


1.5 


1.5 


1.34E+06 


4.15E + 08 


1.26E + 06 


1.11E + 09 


2.29E + 05 


64.87 


-0.13, 1.3 



In Table 2, we summarize the results of simulations where we have varied m^, 
for a fixed value of nih- In all of these cases no jet comes out of the CENBOL (i.e., 
R = 1). In the last column, we listed two spectral slopes a\ (from 10 to lOOkcV) 
and «2 (due to the bulk motion Comptonization). Here, N ms represents the photons 
that have suffered scattering between r g = 3 and the horizon of the black hole. 

In Fig. 5, we show the variation of the spectrum in the three simulations pre- 
sented in Fig. 3(a-c). The dashed, dash-dotted and doubledot-dashed lines are for 
R = 2 (Case la), R = 4 (Case lb) and R = 6 (Case lc) respectively. The solid 
curve gives the spectrum of the injected photons. Since the density, velocity and 
temperature profiles of the pre-shock, sub-Keplerian region and the Keplerian flow 
are the same in all these cases, we find that the difference in the spectrum is mainly 
due to the CENBOL and the jet. In the case of the strongest shock (compression 
ratio R — 6), only 37% of the total injected matter goes out as the jet. At the same 
time, due to the shock, the density of the post-shock region increases by a factor 
of 6. Out of the three cases, the effective density of the matter inside CENBOL 
is the highest and that inside the jet is the lowest in this case. Again, due to the 
shock, the temperature increases inside the CENBOL and hence the spectrum is the 
hardest. Similar effects are seen for moderate shock (R = 4) and to a lesser extent, 
the low strength shock (R = 2) also. When R = 4, the density of the post-shock 
region increases by the factor of 4 while almost 97% of total injected matter (Fig. 
2) goes out as the jet reducing the matter density of the CENBOL significantly. 
From Table 1 we find that the N cen b i is the lowest and Nj et is the highest in this 
case (Case lb). This decreases the up-scattering and increases the down-scattering 
of the photons. This explains why the spectrum is the softest in this case. In the 
case of low strength shock (R = 2), 57% of the inflowing matter goes out as jet, but 
due to the shock the density increases by factor of 2 in the post-shock region. This 
makes the density similar to a non-shock case as far as the density is concerned, but 
with a little higher temperature of the CENBOL due to the shock. So the spectrum 
with the shock would be harder than when the shock is not present. The disk and 
the halo accretion rates used for these cases are rrid = 0.01 and nih = 1- 

In Fig. 6, we show the components of the emerging spectrum for all the three 
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le+07 F 




Fig. 5. Variation of the emerging spectrum for different compression ratios. The solid curve is the 
injected spectrum from the Keplerian disk. The dashed, dash-dotted and doubledot-dashcd lines 
are for R = 2 (Case la), R = 4 (Case lb) and R = 6 (Case lc) respectively. The disk and halo 
accretion rates used for these cases are rri^ = 0.01 and rd^ = 1. See, text for details. 



cases presented in Fig. 5. The solid curve is the intensity of all the photons which 
suffered at least one scattering. The dashed curve corresponds to the photons emerg- 
ing from the CENBOL region and the dash-dotted curve is for the photons coming 
out of the jet region. We find that the spectrum from the jet region is softer than the 
spectrum from the CENBOL. As Nj et increases and N cen b i decreases, the spectrum 
from the jet becomes softer because of two reasons. First, the temperature of the 
jet is lesser than that of the CENBOL, so the photons get lesser amount of energy 
from thermal Comptonization making the spectrum softer. Second, the photons are 
down-scattered by the outflowing jet which eventually make the spectrum softer. 
We note that a larger number of photons are present in the spectrum from the jet 
than the spectrum from the CENBOL, which shows the photons have actually been 
down-scattered. The effect of down-scattering is larger when R = 4. For R = 2 
also there is significant amount of down scattered photons. But this number is very 
small for the case R — 6 as N cen b i is much larger than Nj et so most of the photons 
get up-scattered. The difference between total (solid) and the sum of the other two 
regions gives an idea of the contribution from the sub-Keplerian halo located in the 
pre-shock region. In our choice of geometry (half angles of the disk and the jet), the 
contribution of the pre-shock flow is significant. In general it could be much less. 
This is especially true when the CENBOL is farther out. 

In Fig. 7, the emerging spectra due to the bulk motion Comptonization when 
the halo rate is varied. The solid curve is the injected spectrum (modified black 
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Energy (keV) 



Fig. 6. (a-c): Variation of the components of the emerging spectrum with the shock strength 
(R). The dashed curves correspond to the photons emerging from the CENBOL region and the 
dash-dotted curves are for the photons coming out of the jet region. The solid curve is the spectrum 
for all the photons that have suffered scatterings. See, the text for details. 

body). The dotted, dashed, and dash-dotted curves are for rhh = 0.5, 1 and 1.5 
respectively, rfid = 1-5 for all the cases. The Keplerian disk extends up to 3r g . 
Table 2 summarizes the parameter used and the results of the simulation. As the 
halo rate increases, the density of the CENBOL also increases causing a larger 
number of scattering. From Fig. 4a, we noticed that the bulk velocity variation 
of the electron cloud is the same for all the four cases. Hence, the case where the 
density is maximum, the photons got energized to a very high value due to repeated 
scatterings with that high velocity cold matter. As a result, there is a hump in the 
spectrum around 100 keV energy for all the cases. We find the signature of two 
power-law regions in the higher energy part of the spectrum. The spectral indices 
are given in Table 2. It is to be noted that oli increases with nih and becomes softer 
for high rAh- Our geometry here at the inner edge is conical which is more realistic, 
unlike a sphere (perhaps nonphysically so) in Laurent & Titarchuk (2001). This may 
be the reason why our slope is not the same as in Laurent & Titarchuk (2001) where 
a 2 = 1.9. In Fig. 8, we present the components of the emerging spectra. As in Fig. 
6, solid curves are the spectra of all the photons that have suffered scattering. The 
dashed and dash-dotted curves are the spectra of photons emitted from inside and 
outside of the marginally stable orbit (3r s ) respectively. The photons from inside 
the marginally stable radius are Comptonized by the bulk motion of the converging 
infalling matter and produces the power-law tail whose spectral index is given by 
a 2 (Table 2). 
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Fig. 7. Bulk motion Comptonization spectrum. Solid (Injected), dotted (M^ = 0.5), dashed 
(Mh = 1), dash-dotted (M^ = 1.5). = 1.5 for all the cases. Keplerian disk extends up to 3.1r 9 . 
Table 2 summarizes the parameters used and the simulation results for these cases. 

4. Concluding remarks 

In this paper, we extended the results of our previous work on Monte-Carlo sim- 
ulations (Paper I). We included the outflow in conjunction with the inflow. The 
outflow rate was self-consistently computed from the inflow rate using well-known 
considerations present in the literature (Das et al. 2001 and references therein). We 
compute the effects of the thermal and the bulk motion Comptonization on the soft 
photons emitted from a Keplerian disk around a black hole by the post-shock region 
of a sub-Keplerian flow which surrounds the Keplerian disk. A shock in the inflow 
increases the CENBOL temperature, increases the electron number density and re- 
duces the bulk velocity. Thermal Comptonization and bulk motion Comptonization 
inside the CENBOL increases photon energy. However, the CENBOL also generates 
the outflow of matter which down-scatters the photons to lower energy. We show 
that the thermal Comptonization and the bulk motion Comptonization were possi- 
ble by both the accretion and the outflows. While the converging flow up-scatters 
the radiation, the outflow down-scatters. However, the net effect is not simple. The 
outflow parameters are strongly coupled to the inflow parameters and thus for a 
given inflow and outflow geometry, the strength of the shock can also determine 
whether the net scattering by the jets would be significant or not. Sometimes the 
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Fig. 8. Components of the emerging spectrum for the Cases 2(a-c). Solid curves are the spectra of 
all the photons that have suffered scattering. The dashed and dash-dotted curves are the spectra of 
photons which are emitted from inside and outside of the marginally stable orbit (3r g ) respectively. 
The photons from inside the marginally stable radius are Comptonizcd by the bulk motion of the 
infalling matter. Here the jet is absent. 

spectrum may become very complex with two power-law indices, one from thermal 
and the other from the bulk motion Comptonization. Since the volume of the jet 
may be larger than that of the CENBOL, sometimes the number of scatterings suf- 
fered by softer photons from the electrons in the jet may be high. However, whether 
the CENBOL or the jet emerging from it will dominate in shaping the spectrum 
strongly depends on the geometry of the flow and the strength of the shock. We 
also found that the halo can Comptonize and harden the spectrum even without 
the CENBOL. 
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